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Abstract 

We constrain the mass abundance of nnclustered primordial black holes (PBHs), 
formed with a simple mass distribution and subject to the Hawking evaporation and 
particle absorption from the environment. Since the radiative flux is proportional to the 
numerical density, an upper bound is obtained by comparing the calculated and observed 
diffuse background values, (similarly to the Olbers paradox in which point sources are 
considered) for finite bandwidths. For a significative range of formation redshifts the 
bounds are better than several values obtained by other arguments Q pbh < 1CT 10 ; and 
they apply to PBHs which are evaporating today. 
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1 Introduction 


A variety of observations make a compelling case for the existence of dark matter in 
unknown form(s). Among the discussed possibilities we may mention axions, WIMPS, 
cosmic strings, brown dwarfs and Primordial Black Holes, see JTf|. It has been longely 
recognized that the abundance of evaporating PBHs [|J may be constrained by a variety 
of methods . For example, Carr and Mac Gibbon (see 0 and references therein) used 
the gamma-ray diffuse background to establish constraints. Similarly, Liddle and Green 
0 reviewed a wide variety of methods to study a possible PBH population. Generally 
speaking, the methods dealing with the evaporation of PBHs apply to the final phase of 
these objects only, and therefore do not make full use of the radiative thermal flux at 
larger wavelengths predicted by the Hawking process of quantum evaporation 0,[[|]. 

The formation of PBHs from primordial fluctuations was first studied by Carr J5JJ. 
In that work he established that large, gaussian primordial fluctuations from a scale- 
free primordial spectrum may be responsible for a power-law PBH mass function. A 
new twist to the problem has been added by the recent work of Niemeyer and Jedamzik 
0,0 who used the concepts of critical phenomena to study the initial mass function for 
Gaussian and ’’mexican hat” initial fluctuation spectra. They have argued that in these 
cases PBHs may form with masses below the horizon mass, and found a mass function 
which is not a pure power-law. Black holes may form at all epochs (see Yokoyama 0, 
and Kawasaki and Yanagida [?] for recent work on PBH formation from the collapse of 
large density perturbations in double-inflation and supergravity models) and since there 
is considerable uncertainty about the actual processes that effectively formed PBHs, we 
seek here a general method to constrain the population which may be used independently 
of the specific PBH formation mechanism(s) (see for instance, Hansen et. ah, and Bousso 
0 for a discussion of other possibilities). We shall evaluate a simple (but useful) model 
consisting in a single-mass scale, which may be associated to the typical mass-scale of any 
given mass function showing such a peak (this situation is very well satisfied by models 
in which large perturbations collapse to form PBHs). 

The simplest reexamination of PBH physics has shown that, as a result of particle 


absorption from the expanding background, there is a maximum value of the mass (termed 
’’critical mass” in |TTJ) splitting an evaporating (or subcritical) from a non-evaporating 
(or supercritical) subset of a general mass distribution. The critical mass separating these 
regimes for each redshift is therefore a useful boundary to discuss the physics of the PBHs 
and will be defined in the next section . A detailed quantitative discussion of these issues 


has been given in Refs.|10|,[^]. Due to the absorption of particles from the thermal 
background, PBH evaporation (and therefore some limits to their abundance) have to be 
reexamined. Actually we shall see that the existence of a regime of ~ constant mass in 
the life of PBHs is very important for a full evaluation of their survival. Previous studies 
131 have not considered the regime of quasi-const ant mass and thus need to be updated 


to account for the evaporation at each evaporation scale labeled by the redshift. 


2 Formalism 


In order to simplify our evaluations, let us consider a homogeneous and isotropic cos¬ 
mological distribution of PBHs. The radiative thermal flux from an evaporating black 
hole, as measured by an observer at z = 0 is given by F(M ) = , where Dl(z) = 


H 0 


1 + 


z(l-go) 


y/l+2qoz-\-l+zqo _ 


and go is the deceleration factor (see, for instance, ||). 


The luminosity function L(M) is given by L(M ) = L^Muaw/M) , with L 0 ~ 2.6 x 10 16 erg s” 1 
and it is simply the Stefan-Boltzmann’s law applied to an evaporating black hole. We 
have scaled the mass to a ’’natural” constant Mnaw = 10 15 g (hereafter the Hawking 
mass), since this is the mass of a PBH which is completing its evaporation today. Thus 
we write /x* = and so on along this work. 

The subscript ” * ” will be used for the initial values of any quantity. Moreover, we 
shall restrict our considerations to semiclassical PBHs, i.e. those that satisfy M M p i, 

M p i being the Planck mass. 

We define d/x£(/x, z) as the numerical abundance of PBHs dn at redshift z with masses 
between /x and /j + d/d per unit horizon volume; then £(/x, z) denotes the number of PBHs 
by horizon volume and by mass interval and A Vh or {z)dfi^(n, z) is the number of PBHs 
contained within this volume. An upper limit on the mass density of radiating PBHs 

























is obtained from the requirement SF tota i(PBHs ) < 5Fb ac k, with 5Fb ac k the value of the 
measured background flux in a given frequency interval. 

To begin the calculation we assume a dilute PBH gas, i.e. a negligible collision term 
(PBH-PBH) between these objects. In other words we assume that the cosmic expansion 
rate is much larger than thePBH collision rate. Our aim will be to obtain windows of 
allowed abundances, for which we limit the PBH mass abundance in function of initial 
mass versus formation redshift. 

As discussed in Ref. [[3IJ , the minimal formalism that describes these features is given 
by the set of equations 
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Where g to tai denotes the total density (radiation plus PBHs), Q{t) gives the total heat 
input from the evaporation from the subcritical PBH population (Q{t) = / dmrhN(m,t )), 
R(t) is the usual scale factor, f p bh = fpbhit 1 ! is the distribution function that de¬ 
scribes the cosmological evolution of PBHs, their numerical abundances and higher mo¬ 
menta and f3 denotes the peculiar velocity of the PBH as given by /3 = ( v/c ). In eq.(2) 
A(g) ~ 10 26 g 3 s - 1 h(/i) with /r(/r > 10 2 ) ~ 1 and h(n « 1) ~ 100, describing the degrees 
of freedom of the emitted particles according to the Standard Model. 

The critical mass value (as introduced in |L(|) separating the evaporating (subcritical) 


from the non-evaporating (supercritical) regimes is obtained from the condition /i = 0 
and reads 
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with To ~ 2.7 K the present value of the CMBR temperature, and the other symbols 
have standard meaning. In terms of the Hawking mass we write eq.(5) as fJ> c (T) ~ (t°/t 0 ) • 
Once the individual fluxes from the evaporating subpopulation are calculated, the 
total flux can be obtained by integration 


5F pbh = 


1 r z f dz 


An Jo 


Dl 2 (z) 


A V/i 0r (■*) 

aT - 


rfj’c(z) 

1 6 fipi 


4/i £(ii,z)L(n) 


( 6 ) 


where Zf denotes the formation redshift and 6 is a dimensionless number 0(1). Pri¬ 
mordial Black Holes above n c (z) do not contribute to the radiation since for them the 
absorption term is dominant by definition. As it stands from eq.(6), a suitable evaluation 
of 5F p bh requires integration over the redshift associated to the expansion. This feature 
is depicted in Fig.l which shows the light-cone and the received radiation emitted from 
these objects. All those PBHs with /i* < 1, the Hawking mass, formed at zj 1 had 
evaporated at z evap (ji x ) > 0. 

Generally speaking the functional form of the amplitude £(//, z) will depend on the 
physics of the process that forms the PBHs and can be quite complicated. A reasonable 
ansatz is to factorize the usual volumetric dilution from the mass dependence. Thus, the 
function £(//, z) for the original || Carr mass function is chosen to have the form 


„ N £ C (/b-) ^ h n{Z) 

£(/b-) = 77 77 \ = £o 


V hor {z ) ^V hor (z) i>] 

and, for instance, the critical Niemeyer-Jedamzik is not a power-law and we refer to 

the references 0] , [7] for a full account of their work. A large class of initial mass functions 
are nonetheless included in the ansatz class. 

The single-mass scale mass function to be studied in this paper is simply 


£(v,z) = Co 


- /a*) 
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with the horizon volume given by Vh, or (z ) = dryCr, the numerator denotes the present 

0(^1 + z) 

value for the particle horizon volume and the numerical density in black holes is normalized 















to the present value when z — 0. This normalization factor will be canceled by the same 
factor present in A14 or (^) in eq.(6).The subscript ”c” denotes comoving values and £o is 
the amplitude, interpreted as PBH number by mass interval, i.e. PBHs x g~ l . This work 
will evaluate mass constraints using the mass function of eq.(8), considered as a basic 
model. Other examples like the critical collapse spectrum ||, @ discussed recently, or 
Carr’s power-law spectrum can be worked out, although the evaluations are quite involved. 
However, it is important to stress that all power spectra that produce a substantial number 
of PBHs present sharp peaks, and therefore look pretty much like delta-functions Jl3j . 


To proceed with the evaluation of the bounds we must make a connection between the 
mass function amplitude and the solutions of eqs.(l-4) in the following way: the number 
density of PBHs is given by 


npBH(z) = j dM£(n-> z) (9) 

On the other hand, we may be able to evaluate iipbh(z) using the kinetic formalism 
(see [|1(J) as given by 


9* 

(27T) 3 

with p = mc/3 is the non-relativistic momentum. 

Some algebraic manipulations enable to cast eq.(9) in the form 
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Here, g* is the statistical weight and we integrate up to the critical mass g c (z) since 
above it the evaporation does not take place. Note that for g > g c (z) the third term of 
eq.(4) changes sign and the lower limit 9g p i precludes contribution of quantum PBHs near 
the Planck scale. Now, we may rewrite Hpbh(z) using the same normalization factor as 


rpc(z) 

npBH(z) = M Haw dg£(g,z) (12) 

J6p.pt 

and by comparing eqs.(ll) and (12) we obtain 












£0,z) = 9 * C — J o df3f3 2 f pbh ([3, fi, z) (13) 

Clearly the microphysical framework is a powerful tool, since it can describe all the 
kinematical effects as explicitly displayed in the limits of the integrals eqs.(ll) and (12) 
above. A complete analysis of the solutions is beyond the scope of the present work 
and will be treated elsewhere. Here, we will analyze only the radiation emitted by these 
objects. 

Our general picture is now formally complete, in the sense that solving the formalism 
f° r fpbh(/3, H, z) we would be able to set upper limits to the abundance of PBHs using 
the thermal emission of any mass spectrum. Since the work done in the literature mainly 
addresses the final explosive phases of these objects, which contributes to the gamma-ray 
band and to the cosmic ray flux,(see for example J3j] and references therein), we may view 
this work as an attempt to generalize those results. 

3 Application to a single-scale mass function and 
comparison with other methods. 

To proceed we must now evaluate the last integral found in eq.(6), for the delta mass- 
function case where £(/i, z ) oc — /i*). The limits on £ 0 will thus apply to an unclustered 
isotropic distribution of PBHs that are subdominant for the total mass balance of the 
universe. These hypothesis are in fact not very restrictive since for low peculiar velocities 
and the masses considered, the time scale for clustering is much larger than the Hubble 
time Hq 1 . 

As explained above, we will assume that at Zf the mass is concentrated at one discrete 
value. Figure 2 shows the two important regimes to be considered here: /x* > 1 or //* < 1 
in terms of the formation redshift and the critical mass fj, c {z). To visualize the differences 
we shall discuss them separately. 

First regime: /i* > 1 

In this case, all PBHs formed at Zf are still evaporating today. Here, M Haw 10 15 g is 
the Hawking mass, i.e. it is the mass that completely evaporates today. We recall that the 



time scale for evaporation depends only of //* and it is given by t evap = /(ffo)^o _1 (f t *) 3 - 
We may split this regime in two sub-cases depending on wether Zf the PBH mass is larger 
or smaller than the critical mass (n(zf) > Hc( z f) or fi(zf) < fi c (zf)) at Zf. In the first 
subcase, the PBHs the onset of evaporation will be delayed until they cross the critical 
mass curve, as will be described later. Now, we will study the subcase fi*(zf) < /i c {zf) 
which consists of initially subcritical PBHs. 


3.1 Subcase /i*(£/) < /i c (zj) 

These PBHs will start to evaporate immediately after their formation. The evolution of 
this spectrum with z will be given by 


£{n{z),z< z f ) = £ 0 


W) ~ MAP 
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where n(z) is given by the /i{z) = /x*F(/x*, z ini , z), and the evaporation function 
F(/i*,z ini ,z) is 


Zi n i , z) ~ 



(15) 


which can be approximated by (1 — — 3 ^ + ^ 3 / 2 )^ 3 ~ 1, for z ini >> 1. 

This form is valid as long as go = 1/2, because we have used that the horizon volume 
at z is related to the horizon volume at Zf by Vh or (z) = Vhor( z f))- We do not 

consider other cases for g 0 in this article, but we checked that the dependence of the 
results with g 0 is very mild. 

Substituting these relations into eq.(6), we evaluate the radiation received today (z = 
0) from all these PBHs yielding 


= 3 LqMh'J* p z)]~ 

P 4vr 4 D L 2 (z)(l + z) 1 ' 

Assuming that PBHs have black body spectra, we can apply Wien’s law to the mass- 
temperature relation of black holes: T(/i) ~ 2 x 10 n /i/(/i) in order to obtain an asso¬ 
ciated wavelength of the maximum A(/i*)T(/i*) = 0.29 cm K. Since we are considering 



the radiation emitted by PBHs at cosmological distances, we must correct for the ex¬ 
pansion of the received wavelength. We know that these relation would be satisfied by 
A obs(z) = A Substituting into the mass-temperature relation 

A obs{z = 0, n(z')) = A + z') (17) 

And considering that at any time the mass is given by /x(T) = /x*F(/x*, z ini , z') we 
obtain 


Ad; /X*) A 0 /X*F(/X*, Zi n i, a, ) (1 -|- z ) (18) 

The prefactor in front of (1 + z') is the emitted wavelength when the initial PBH with 
mass /i* had the mass /i{z) at the correspondent temperature, and A 0 ~ 1.5 x 10~ 12 cm is 
a reference scale. Then, we recover z ini , z') = By its very definition 

F(/x*, z ini , z’ = z ini ) = 1, and thus 
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After all these manipulations the total flux obtained is given by 
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where z ini = z f if /x* < fJ, c (z f ) or z ini = z cross (n*) if /x* > /x c (z/). 


An inspection of eq.(20) shows that there is a potential mathematical divergence at 
the lower limit of the integral. Therefore we have integrated formally the expression up 
to a minimum parameter e Q , which acts as a cutoff. The integral I(e e ,Zf) = f** dz ) 
admits the expansion 


I(e Q ,z f )~Y( qo )(H 0 /c) 2 
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with V(qo) ~ 0(1) and e© << 1. The question arises about how this e© affects the 
final results. We shall show below that a careful consideration of the different regimes of 



the PBHs makes unnecessary the imposition of any e Q , and in practice there is no physical 
divergence of the quantities, contrary to the naive expectation from eq.(20). 

Now, if we impose that 5F pbh < 8F back , we obtain a constraint on the amplitude £ 0 
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Although the actual limits may be considered for any small interval, we have defined 
in practice an average background value for the radiation over a large portion of the 
spectrum as 
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.10 6 ergs~ 1 cm~ 2 . 

Some care has yet to be taken to take into account the full dependence of the mass 
density in PBHs g pb h(z) inside the horizon volume. After some algebric manipulation and 
using the properties of the delta function we derive the expression 


_ f 1 + ^ \ 3 €o(M Ha w) 2 w N 

Q pbh (z)- { 1 + Zf ) V hor (z f ) (/ ^^* 24 

which is valid for any z > z evap (n *). If we divide eq.(24) by g c (z) = £> c (0)(l + z) 3 
(again for q 0 = 1/2), we will obtain Q pbb (z) for any z which satisfies the causal relation 
(that is, fi pb h A giior) i 


(1 + Zf) < 


4.24 x 10 28 



(25) 


This expression precludes the existence of PBHs with masses larger than the one 
contained within the cosmic horizon, i.e. g pb h(z) < ghor(z), but admits the case g, pb h -C 
Hhor as advocated in Refs. [6,7]. For n phh ~ 10 11 , the critical mass value today, we have an 
upper limit at which PBHs can be formed still satisfying eq.(25), namely Zf ~ 10 16 . 
To estimate actual bounds we evaluate the expression at z — 0. Then, we have 
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Now, using eq.(23) and dividing it by the critical density g c (z = 0) ~ 2 x 10 29 gem 3 , 
we have (imposing Zi n i = Zf) 


tt p bh{ 0) < tt p bh(z f , n*) = 8.7 x 10 8 (1 + z f ) 2 n* 3 F(n*, z f , z = 0 )e o < SF back > (27) 

where F(/r*, Zf, z = 0) rs_/ 0 ( 1 ). 

The interpretation of eq.(27) is the following: the right hand side sets the maximal 
abundance in PBHs today ((z = 0)) allowed by the sky brightness assumed to be fully 
produced by the evaporation of PBHs formed at Zf with initial masses given by /r*. These 
PBHs were born above the Hawking mass and therefore they are evaporating today. We 
have also required that their masses were initially below the critical mass value today 
~ 10 26 g, otherwise our constraints do not apply in this form (see below). From now, we 
shall denote the product e Q < SF ba ck > which appears conspicuously by 6b- 

Our bounds based on the sky brightness will be tighter (by construction) inside a finite 
region in the parameter space defined by /j* versus Zf. This finite region corresponding to 
ft P bh(zf, //*) < 1 CT 10 may be further explored to understand which is the physics allowing 
these bounds. Taking the logarithm on both sides of eq.(27) we find 
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with es = e© < dFback > as stated above. This inequality displays a relation between 
the initial mass /i* and formation redshift Zf bounding the ” coolest PBH population”, 
since any PBH population below this line will be contributing so much to the sky bright¬ 
ness that we can certainly rule it out . This curve is a consequence of the Doppler effect, 
since for Zf above this curve, the cosmic expansion diminished the total flux below the 
background value and big masses correspond to very cold PBHs, for which the constraint 
is necessarily very poor. It is clear that the ^/-threshold is mass-dependent and therefore, 
the curve must drop at later times. 

In terms of the mass we may write eq.(28) in the form 
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Meaning that all those PBHs formed at Zf satisfying both n* > 1 and the eq.(28) 
above are much more scarce than Q p bh(n*) < 10~ 10 , for a given combination value of the 
Zf and the parameter e B . 

The very fact that 1 < /r* < leads to some additional constraints avoiding the 

imposition of an e©. Actually, since fx c {zf) = /z c (0)(l + z/) -1 , we note that these relations 
can hold only for (1 + Zf) < 10 11 . But the condition /i* > 1, eq.(28) requires also that 
e©< SFback > < 10 _3 (1 + ZfY' 2 . Taking the logarithm as before 

-\loge B > 1 + \log(l + z f ) (30) 

Using that (1 + Zf) < 10 11 we obtain an upper bound to the value of the parameter 
<7b 


e B < 10 - 25 (31) 

Eq. (31) is a direct consequence of imposing Q p bh(l-i*, Zf) < 1CT 10 for //* < 1 and 
Zf 1. We see that e B is bounded from above by a very small number, in fact so small 
that its associated distance D B (e Q ) ~ 1 cm. This means that we can exclude PBHs in the 
discussed conditions regardless of the irrelevant value of e Q , which is a ”physical” zero. 

Analogously, we must satisfy /r* < fj, c (z = 0), which implies e B > 10~ 59 . Fig.3 
displays the region between logzf and loge B in which all those PBHs satisfying 1 < /x* < 
(i + 7 )°' 3 (€b) i /3 are more scarce than ^^(^u*) < 10 ~ 10 . 

The region in the parameter space logzf versus logg * where Q p bh(g*) < 10 _1 ° is then 
bounded by the following inequalities 

2log(l + Zf) + logn* < 57.2 (32) 

log fi * < logfj, c ( 0) ~ 11 (33) 


logzf > logz min ~ 4 


(34) 



log/d * > 0 (35) 

and eq.(28). It is clear that eq.(32) describes the causality constraint, i.e. there can 
be no PBH with mass greater than the horizon mass at zj. The eq.(33) states that our 
method applies to evaporating PBHs only. Eq.(34) says that PBHs do not form as late 
as z < 10 4 , since we know that the environment must be very smooth at late times (this 
arbitrary prescription can be relaxed straightforwardly). Finally, eq.(35) describes the 
choice /i* > 1. The compact region where the abundance of PBHs is < 10” 10 in this 
plane is displayed in Fig.4. 

3.2 Subcase /i c (^/) < /i* < Hhor( z f) 

The range of masses implies (1 + Zf) < 4.24 x 10 28 , therefore we can rewrite the previous 
expressions to obtain for £ 0 

£o < 6.8 x 10 52 (1 + z cross (id*)) 2 id* 2 h 0 ~ 2 e B (36) 

And considering that (1 + z cross fi*) 2 ~ 10 22 /i*~ 2 we insert these relations into eqs.(24) 
and (31) to yield 

0, f-i* ^ Rc(^/)) ^ 8.7 x 10 i-i^eqF (/x*, z 0, z^^) (37) 

Requiring that < 10~ 10 as before, the corresponding masses must satisfy 

10~ 25 

/J'c(zf) < /i* <- (38) 

Pb 

Analogously to the former case /i* > 1 then the cutoff is limited to 

< 10- 25 , (39) 

And the same considerations as before apply. Since /x* < g c (z = 0) ~ 10 11 must also 
hold, we obtain the limit 


e B > 10 ” 34 , 


( 40 ) 



again corresponding to microscopic lengthscales, and therefore not relevant for the 
bounds so obtained. Finally, requiring that 1 °^ > /j, c (zf, we obtain the inequality 

^ - 36 <«) 

Note that eqs.(38) and (39) are completely independent of Zi n i because the latter 
dependence cancels out. Then, it is eq.(39) the one that sets the maximum value for the 
lower cutoff in redshift, already shown to be irrelevant. Although we may have considered 
the subcase /i* > gfzf in two parts, (the first 1 > g c (zf and the second 1 < gfzf), 
they actually differ only in the allowed value of zf. for the first we have (1 + zf > 10 11 
and for the second case, (1 + zf) < 10 11 . 

Figure 5 displays a region in the plane log( 1 + zf versus e# in which we find PBHs 
satisfying Q pbh (gf < 10~ 10 and for which yU* > gfzf. All these PBHs satisfy the 
conditions yU* > 1, /x* < yU c (0) ~ 10 11 and logzf < 28.6 as well (for logzf > 28.6 we would 
have g hor < 1). 


Second Regime: yU* < 1 


When we deal with PBHs born with masses smaller than the Hawking mass initially, 
we may substitute the lower limit of the integral eq.(21) by z evap (gf, because no PBH 
will survive its own timescale for evaporation. Therefore the integral I(z evap (gf , z ini ) = 
fzTa P (n*) dz D*0) is explicitly evaluated as 
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where G(gf = 
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Following the same reasoning as in the previous case we obtain for £ 0 
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where z mt falls in one of these two subcases: when yU* > g c (zf we have z m i 
and when yU* < g c (zf we have z„ a = Zf. 


(43) 
10 12 /n. 



Using the expressions above we finally find, after substituting eq.(43) into eq.(31) and 
dividing by Q c (z) 
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Since that F(fi^, z ini , z = z evap (ix *)) = 0, the result eq.(44) describes correctly the 
evolution of the abundance from Zf until z evap (^*), taking into account the Hawking 
evaporation. The most interesting case obtained from eq.(44) applies to z = Zf, from 
which we derive the maximum abundance allowed at formation still consistent with the 
sky brightness today < 5F back >) 
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for /i* < and 


^pbh(Zf) < 9 x 10 
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_G(/i*) + /n(10 n p* 

for p* > Hc(z f ). 

However, when we try to force eqs.(45) and (46) to yield values Q pb h < 10 ~ 10 we did not 
find physical windows for reasonable values of PBH masses. The reason is that to obtain 
Qpbh < 10 “ 10 ; p* have to be very small and thus violate our semi-classical assumptions. 
This limitation is to be expected, since those PBHs evaporated a long time ago when the 
cosmic radiation was ~ (1 + z) 4 times more intense than the present background. This 
means that we can not use the cosmic isotropic radiation in order to estimate constraints 
for the sub-Hawking PBH population, and other methods (as described, for example, by 
Green et al. 0) must be used to obtain useful limits to their abundance. 


4 Conclusions 


We have shown that under a set of circumstances a limit on Vt pbb better than 10 10 can 
be obtained by using the background brightness of the sky. This is quite stringent if we 












compare it with the values given by other methods |4j]. We attempted to evaluate the 
evolution of maximal abundance Q p bh(z) for any 0 from observations performed today 
(z — 0), a process that requires integration over z (and eventually over the actual mass 
spectrum from whatever the formation process, not attempted here). 

The good news here is that PBHs can be excluded in a fairly large window (as seen 
in Fig. 4) if they were born above the Hawking mass (but below 10 26 g). We see four 
physical reasons for this fact: 

1) Many of these PBHs still evaporate today, and therefore contribute to the back¬ 
ground in various bands. 

2) The cosmic radiation today is much less intense than at high redshift values, and 
therefore constitutes a useful natural background. 

3) The cosmic expansion today is quite weak, then the Doppler damping of the emitted 
radiation does not have any killing effect. 

4) PBHs initially above the critical mass at some Zf delay the start of their evaporation, 
therefore they injected energy later, when the expansion is less effective to damp the 
radiation. 

The constraints so obtained are weakly dependent of g 0 or the cosmological constant, 
and the mass range probed by this method is larger than previous works. Even though 
we take into account the existence of the critical mass and the redshift of the emitted 
radiation through the function I(e & ,Zf), we found that our analysis renders useful limits 
for PBH masses initially above ~ 10 15 g, but not for those below it. In addition, PBHs 
above 10 26 g must be limited by other methods because they do not evaporate at all. 
The consideration of the different cases leaded us to conclude that e 0 need not to be 
imposed by hand. Moreover, it happens to be an extremely small number in all cases, 
with microscopic associated distances D^(e 0 ). Physically, this means that the spatial 
distribution of PBHs does not need to satisfy any specific requirement near the earth, 
and thus the bounds so obtained are quite general. 

The upper limits to the number density of PBHs were obtained from the simple re¬ 
quirement 6F p bh < SFb ac k, since other astrophysical mechanisms may be also contributing 
to the cosmic background brightness as measured by our detectors. The limits have been 


derived for a Dirac’s delta mass function, whose main features are much simpler than any 
other choice. We can think this case as a first approach to the general problem, likely 
accurate for peaky distributions of a more general type. Detailed calculations of the latter 
remain an interesting problem for future investigations. 
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Figure captions 


Figure 1. The fate of primordial black holes in the redshift vs. distance-luminosity plane. 
A PBH formed at Zf orm evaporates and the emitted radiation arrives at the observer when 
its path ( explicitely indicated) crosses the vertical axis. If the PBH mass is initially 
below the Hawking mass, then , the corresponding redshift for complete evaporation z evap 
is located at z > 0. 

Figure 2. Mass scales vs. redshift of formation plane. The critical mass curve and the 
Hawking mass constant are explicitely indicated. Note that both cross at z rsj 10 11 . In 
the region A hot (evaporating) PBHs begun to evaporate right after their formation. In 
the region B the PBHs did not begin to evaporate at formation, but only after crossing 
the critical mass curve and entering into the region A. In the region C PBHs were formed 
at very high redshift, greater than z = 10 11 and followed a similar history. Below 
Mjjaw ~ 10 15 g, PBHs that evaporated completely today. Our method of analysis and 
bounds thus obtained are useful for the PBHs formed at the regions A, B and C (termed 
’’First Regime” in the text) but not below the Hawking mass (see text). 

Figure 3. The region inside which i'lpbh < 1CT 10 in the es vs. Zf (hatched) for those 
PBHs born in region A. The boundaries of this region are defined in the text. 

Figure 4. The compact region in which the Vt v bh < 1CF 10 in the redshift vs. mass plane 
(see text). At logg* ~ 11 the PBHs have a mass of 10 26 g, the critical mass value today. 
The inclined straight line bounding region C from above is set by the Doppler effect of 
the cosmological expansion. We have also drawn the horizon mass line and explicitly 
indicated the non-causal region above it. The regions A, B and C correspond to the 
regimes and subcases defined in the text. 

Figure 5. The same as Fig. 3 for the regimes corresponding to the regions B and C. The 
region B spans redshifts from 0 up to 10 11 , when the critical mass equals the Hawking 
mass. The region C corresponds to any redshift value bigger than 10 11 up to 10 28 , the 
moment when the horizon mass equals the Hawking mass (of academic interest only). 



The upper diagonal line reflects the combination of the Doppler effect and the value of 
€b (see text). 
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